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Abstract 

We report the first observation of charmed mesons with the ZEUS detector at HERA 
using the decay channel —> K“7r'’')7r^ (+ c.c.)- Clear signals in the mass 

difference AM=M(D*)-M(D°) as well as in the M{K'k) distribution at the mass are 
found. The ep cross section for inclusive production with < A GeV^ in the yp 
centre-of-mass energy range 115 < lU < 275 GeV has been determined to be (32±777) nb 
in the kinematic region > 1.7 GeV, |7/(D*)| < 1.5}. Extrapolating outside this 

region, assuming a mass of the charm quark of 1.5 GeV, we estimate the ep charm cross 
section to be a{ep —> ccX) = (0.45 ± 0.1l7o22) at ^/s = 296 GeV and (W) = 198 
GeV. The average yp charm cross section a{'yp —> ccX) is found to be (6.3 ± 2.2]'^3;g) pb 
at {W) = 163 GeV and (16.9 ± 5.2lg^5®) ph at (IT) = 243 GeV. The increase of the 
total charm photoproduction cross section by one order of magnitude with respect to low 
energy data experiments is well described by QCD NLO calculations using singular gluon 
distributions in the proton. 
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1 Introduction 


In high energy ep collisions at HERA the leading order QCD contribution to charm production 
is the photon-gluon fusion (PGF) mechanism ii- In this process the photon participates 
as a point-like particle {direct photon process) interacting with a gluon from the proton and 
giving a cc pair {'-jg —> cc). The PGF cross section for ep —> ccX behaves like da/dQ'^ ~ Q~^ 
and is dominated by the exchange of almost real photons ~ 0), i.e. by photoproduction 
events where the electron is scattered by a small angle. As a consequence, the contribution of 
the Deep Inelastic Scattering (DIS) regime, ~ 4 GeV^, where the scattered electron is seen 
in the main detector, is expected to be small compared to photoproduction. 


Apart from the direct channel, charm photoproduction at HERA can also proceed via the 
resolved photon processes ii. where the photon behaves as a source of partons which can 
scatter off the partons in the proton (mainly gg —> cc). QGD-based models predict that these 
types of processes dominate over the direct processes for light quark production i, i. This 
prediction has been confirmed by measurements of the ZEUS §0 and HI collaborations. 
The predicted open charm cross section at HERA has two major uncertainties: the mass of the 
charm quark {me) @.0 and the structure functions of the proton and the photon. Next to 
leading order (NLO) corrections have been calculated and found to be substantial |^, |^. The 
full NLO cross section (j{ep —>• ccA) at HERA, obtained using me = 1.5 GeV and the structure 
function parametrisation MRSD'_ |]^ for the proton and GRV HO [|^ for the photon is 
predicted to be ~ 0.9 ph at a/s = 296 GeV. A variation of me by ±0.3 GeV changes 
the values of the cross sections by a factor of 2. The predicted relative amount of the direct 
and resolved contributions depends on the photon structure function parametrisation. For 
DG [M, GRV or AGFGP 1^ the resolved contribution is less than 30%, but if the LAGl K 


parametrisation is used the resolved contribution can be 50% or more and the predicted cross 
section can increase by almost a factor of 2. As a consequence, estimates of the total charm 
cross section a{ep —> ccX) at HERA vary between 0.3 ph and 2 ph. 


We search for open charm production at HERA with the ZEUS detector by looking for the 
fragmentation products of the heavy quarks which produce a D*(2010)^. The method relies on 
the tight kinematic constraints of the decay chainQ 


D*+ 




(K 7r^)7rc 


( 1 ) 


where the momentum of the pion coming from the D 
D* rest frame. Gonsequently, the mass difference |]I7 


(‘soft pion’, 7is) is just 40 MeV in the 


AM = M(D°7rs) - M(D°) = 145.42 MeV 


can be measured much more accurately than the D* mass itself. This low Q-value of the 
D* —> D°7r5 decay yields a prominent signal in an otherwise phase space suppressed kinematic 
region, the threshold of the M(K7r7r5') — M(K7r) distribution. 

We assume that the fraction of D*^ originating from bb is negligible [Q. 

^In this analysis the charge conjugated decay chain D*“ —> D^tt^ —> (K+7r“)7rg is also included. 
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2 Experimental Setup 


Data were collected during the 1993 running period, when protons of energy Ep = 820 GeV 
were colliding with electrons of = 26.7 GeV. Gollisions took place between 84 electron and 
proton bunches with typical beam currents of 10 mA. Additional unpaired bunches of electrons 
and protons allowed an estimation of beam related background. 

The total 1993 luminosity collected by ZEUS was ~ 550 nb“^, of which 486 nb“^ were used 
in the present work. This restricted sample contains runs taken with stable trigger conditions 
and the tracking chambers operating fully in the nominal magnetic held. 


Gharged particles are measured by the ZEUS inner tracking detectors, which operate in a 
magnetic held of 1.43 T provided by a thin superconducting coil. Immediately surrounding 
the beampipe is the vertex detector (VXD) |18] consisting of 120 radial cells, each with 12 
sense wires. It uses a slow drift velocity gas and the presently achieved resolution is 50 jum in 
the central region of a cell and 150 pm near the edges. Surrounding the VXD is the central 
tracking detector (GTD) which consists of 72 cylindrical drift chamber layers, organised into 9 
‘superlayers’ [^. With our present calibration of the chamber, the resolution of the GTD is 
around 260 pm. The resolution in transverse momentum for tracks going through all superlayers 


is ~ \J (0.005)^pf^ + (0.016)^ where pt is in GeV. The single hit efficiency is greater 

than 95%. The efficiency for assigning hits to tracks depends on several factors: very low px 
tracks suffer large systematic effects which reduce the probability of hits being attached to them, 
and the 45 degree inclination of the drift cells also introduces an asymmetry between positive 
and negative tracks. Nevertheless, the track reconstruction efficiency for tracks with > 0.1 
GeV is greater than 95%. Using the combined data from the VXD and GTD, resolutions 
of 0.4 cm in Z and 0.1 cm in radius in the XY planeQ are obtained for the primary vertex 
reconstruction. From Gaussian hts to the Z vertex distribution, the r.m.s. spread is found to 
be 10.5 cm, in agreement with the expectation from the HERA proton bunch length. 


The high resolution uranium-scintillator calorimeter (GAL) is used in the present ana¬ 
lysis to calculate global quantities of the events. It covers the polar angle range between 
2.2° < 9 < 176.5°, where 0 = 0° is the proton beam direction. It consists of three parts: the 
rear calorimeter (RCAL), covering the backward pseudorapidityQ range (—3.4 < rj < —0.75); 
the barrel calorimeter (BCAL) covering the central region (—0.75 < p < 1.1); and the forward 
calorimeter (FGAL) covering the forward region (1.1 < p < 3.8). The calorimeter parts are 
subdivided into towers which in turn are subdivided longitudinally into electromagnetic (EMG) 
and hadronic (HAG) sections. The sections are subdivided into cells, each of which is viewed 
by two photomultiplier tubes. Under test beam conditions the GAL has an energy resolution, 
in units of GeV, of ax = 0.35^jE{GeV) for hadrons and ax = 0.18Y^i7(GeV) for electrons. 
The GAL also provides a time resolution of better than 1 ns for energy deposits greater than 
4.5 GeV, and this timing is used for background rejection. 


We use two lead-scintillator calorimeters (LUMI) pT[| to measure the luminosity as well as 
to tag very small processes. Bremsstrahlung photons emerging from the electron-proton 
interaction point (IP) at angles below 0.5 mrad with respect to the electron beam axis hit 

^The ZEUS coordinate system is defined as right handed with the Z axis pointing in the proton beam 
direction, hereafter referred to as forward, and the X axis horizontal, pointing towards the centre of HERA. 

^The pseudorapidity rj is defined as — ln(tan |), where the polar angle 9 is taken with respect to the proton 
beam direction. 
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the photon calorimeter 107 m from the IP. Electrons emitted from the IP at scattering angles 
TT — 0' < 6 mrad and with energies 0.2 E^. < E'^ < 0.9 E^. are deflected by beam magnets and 
hit the electron calorimeter placed 35 m from the IP. 


3 Trigger Conditions 


Data are collected with a three level trigger The hardwired First Level Trigger (FLT) 

is built as a deadtime free pipeline. The FLT decision is based on energy deposits in the 
calorimeter and luminosity detectors, on tracking information and on the muon chambers. We 
require a logical OR of hve conditions on sums of energy in the calorimeter cells: either the 
BCAL EMC energy exceeds 3.4 GeV; or it exceeds 2.0 GeV, if any track is found in the GTD; or 
the RGAL EMG energy, excluding the towers immediately adjacent to the beam pipe, exceeds 
2.0 GeV; or it exceeds 0.5 GeV, if any track is found in the GTD; or the RGAL EMG energy, 
including the beam pipe towers, exceeds 3.75 GeV. 


The Second Level Trigger (SLT) uses information from a subset of detector components 
to separate physics events from backgrounds. It rejects proton beam-gas events using particle 
arrival times measured in the forward and rear sections of the calorimeter, reducing the FLT 
trigger rate by almost an order of magnitude. 

The Third Level Trigger (TLT) has available most of the event information. It is used to 
apply stricter cuts on the event times and to reject beam-halo and cosmic muons. Beam-gas 
interactions are rejected by requiring: 


• a reconstructed Z vertex position within 75 cm of the nominal interaction point, 

• -pzi) > 3 GeV, 

• EiPzi/EiEi <0.9 , 

where the sums run over all calorimeter cells i and pzi is the Z-component of the momentum 
vector assigned to each cell of energy Fj. 

The following additional TLT requirements are made in order to further reduce the back¬ 
ground and the output rate: 


• iPT)max > 0.5 GeV and {pT)max > 0.5 GeV and (py)max + {PT)max > 1-3 GeV, where 
{PT)max is the track of positive or negative charge with the highest pt, 

• transverse energy outside a cone of 10° around the proton direction in excess of 5 GeV; 
or J2i{.Ei — pzi) > 15 GeV; or Y^iPzilHi Ei < 0.8; or an electron with energy larger than 
5 GeV detected in the LUMI electron calorimeter. 
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4 Kinematics 


Neutral Current ep scattering e {k) + p (P) —^ e {k') + X can be described with the following 
kinematic variables: 


/s = ^{k + Pf ^ ^4PpPe = 296 GeV, 
the total ep centre of mass energy; 

= {k - k'f , 

the four-momentum transfer squared carried by the virtual photon; 

q ■ P 

^ ^ k-P' 

the Bjorken variable describing the energy transfer to the hadronic system; and 

W'^ = {q + PY = -Q^ + ys + ml, 

the centre of mass energy squared of the yp system, where rup is the mass of the proton. 

The variable y can be determined to a good approximation from the hadronic system using 


the Jacquet-Blondel expression [23 


yjB — 


Ei{Ei-pzi) 

2 ■ E, 


with the sum running over all calorimeter cells i except for those associated with the scattered 
electron. 

If the scattered electron is seen in the main detector (DIS events) or in the LUMI electron 
calorimeter (tagged photoproduction events), the variable y can also be obtained from: 

^ P' 1-COS0' 

Ve = 1-^ 


5 Monte Carlo Simulation 


The Monte Carlo programs HERWIG and PYTHIA are used to model the hadronic final 
states in cc production. Both generators include parton showers in the initial and final states. 
Fragmentation into hadrons is simulated with the LUND string model as implemented in 
JETSETp^ in the case of PYTHIA, and with a cluster algorithm in the case of HERWIG. 
The lepton-photon vertex is modelled according to the Weizsacker-Williams Approximation 
^ in the case of PYTHIA, whereas HERWIG uses exact matrix elements for the 


(WWA) 

photon-gluon fusion (PGF) component and the Equivalent Photon Approximation (EPA) 
for resolved processes. 

For these Monte Carlo models a large number of cc events was generated, but only those 
containing at least one charged D*, decaying into D^tts with subsequent decay D° —>■ Ktt, were 
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processed through the standard ZEUS detector and trigger simulation programs and through 
the event reconstruction package. With both Monte Carlo programs we have generated events 
with both direct and resolved components, setting = 1.5 GeV. The parton densities in the 
proton were described by MRSD'_ and by GRV HO in the case of the photon. For systematic 
checks, we also generated events using different parametrisations for the proton (MRSDq ||12|| , 
GTEQ2M' 13 and GRV HO 0) and the photon (DG, AGFGP and LAGl). To check the 


dependence of the results on the charm mass assumed in the Monte Garlo simulations the whole 
analysis was repeated using the default values for rric in both Monte Garlo generators (1.35 GeV 
in PYTHIA and 1.8 GeV in HERWIG). 


6 D* Observation 


6.1 Reconstruction Method 

In order to select a kinematic region where the efficiency of the tracking detectors is high 
and systematic uncertainties are well understood, the following requirements on the tracks are 
made: 


• Pt > 0.16 GeV; 

• |? 7 | < 1.5, corresponding to 25° < 6 < 155°. 


Pairs of these tracks with opposite charge are combined and considered in turn to be a kaon 
or a pion. The combination is accepted as a possible D° candidate if the Kvr invariant mass lies 
in the range: 

1.80 < < 1.93 GeV 


(the nominal value of the D° mass is 1.865 GeV |^). To reconstruct D* mesons, these D*’ 
candidates are combined with an additional track having opposite charge to that of the kaon. 
Assuming this third track to be the soft pion, the mass difference AM = M(K7r7r5) — M(K7r) 
is then evaluated. 


Monte Garlo studies show that after these cuts more than 95% of the decay products of the 
D° satisfy pr > 0.5 GeV, with a mean value of 1.5 GeV. These high transverse momentum 
tracks have a higher reconstruction efficiency and a better track extrapolation to the vertex. 
As a consequence, the following more stringent cut is applied to them: 


Pt{K),pt{7i) > 0.5 GeV. 


From Monte Garlo studies we hnd that the tts travels essentially in the same direction as 
the D* itself. Therefore the \p\ < 1.5 cut on the single tracks limits the //(D*) range to: 

\v{B*)\ < 1.5. 

This cut was thus also applied. 

Moreover, we have restricted our analysis to: 

Pt(D*) > 1.7 GeV, 

since more than 95% of the D*’s fulhl this condition after the above cuts. 
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6.2 AAf and Af(D°) Signals 


DIS events are defined to be those having an electron identified in the CAL with ye < 0.7. The 
for these events is larger than 4 GeV^. We hnd that 20% of the events in which we hnd a D* 
candidate fulhl this requirement. This relatively large fraction of DIS candidates, reproduced 
by the Monte Carlo programs, is due to the higher trigger acceptance for DIS events than for 
photoproduction events, and to the harder pt(D*) and Pt{t^s) spectrum for these events. We 
show the AM distribution for the DIS candidates in Fig. ^ A clear D*^ signal around the 
nominal value of AM is observed which is evidence for D*^ production at HERA in DIS with 
> 4GeV^ 


The 80% of the D* candidates which are not identihed as DIS events have < A GeV^ and 
are called photoproduction events. Of these, 27% are tagged in the LUMI electron calorimeter, 
in agreement with Monte Carlo simulations and the ZEUS photoproduction measurements |^. 
To reduce possible background from DIS, where the electron is not identihed, we require 
yjB < 0.7 1^. The true W and y are underestimated using the Jacquet-BIondel method. We 
correct the measured Wjb with Monte Carlo methods [^, resulting in a corrected W value 
which will be used in the following. Comparison with the W measured from events with a 
LUMI tag shows that the estimated true W has a systematic uncertainty of less than 10%. 
The cut on yjs corresponds to IF < 275 GeV. Furthermore we will restrict ourselves to 
lU > 115 GeV where the acceptance is above 8%. The AM distribution for these photopro¬ 
duction events obtained with the set of cuts described above is shown in Fig. ^ja. A clear peak 
around AM = 145.5 MeV is observed. 


In order to check the background shape, pairs of tracks with the same charge are selected 
for calculating the Ktt invariant mass (wrong charge combination method). This distribution 
is htted with the maximum-likelihood method using the function: A x (AM — where A 

and B are the free parameters of the ht. 

The signal distribution is then htted assuming this function for the background plus a Gaus¬ 
sian to parametrise the signal shape. The corresponding hts and the normalised background 
are also shown in Fig. |^. The background shape parameters obtained in the signal ht agree 
with the values obtained by htting the wrong charge distribution. 

We observe 48 ± 11 D*’s above background, with a signal to background ratio of about 1.5. 
The mean value of the AM signal obtained from the ht is AM = (145.4 ±0.2) MeV, consistent 
with the nominal value. The corresponding width is (0.9T0.2) MeV, in agreement with 0.7 MeV 
obtained from Monte Garlo simulation. The mean W for these events is {W) = 198 GeV. 

To check whether the D*’s are produced according to the decay channel (|ID, we show in 
Fig. |b the M(K7r) distribution for the events in the AM range from 142 MeV to 149 MeV. 
A clear signal is seen around the nominal value of the D° mass. In order to ht this distribution 
we have used an exponential background shape plus a Gaussian for the signal. The number of 
observed D'^’s is 43 ± 12, consistent with the number of D*’s obtained from the ht to the AM 
distribution. We obtain a mean value of M(D°) = (1.854 ± 0.005) GeV, slightly below the 
nominal value, with a width of (18 ± 4) MeV, consistent with the value of 17 MeV obtained 
from Monte Garlo simulation. 


6 
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Cross Sections 


7.1 ep Cross Section 


The cross section is obtained using the expression: 

N 

C X BR X Acc' 

where iV = 48 is the number of observed D*’s, C = 486 nb~^ is the integrated luminosity, 
and BR = (2.73 ± 0.11)% is the combined branching ratio of the decay chain (|l]). The 
acceptance Acc was calculated as the number of detected over generated D* —>• Ktitis decays 
in the kinematic range chosen, using Monte Carlo methods including the trigger simulation. 
We have used the PYTHIA Monte Carlo prediction with MRSD'_/GRV HO structure function 
parametrisations for the proton/photon to correct our data and quote cross sections. The 
overall acceptance in the kinematic region {p 7 ’(D*) > 1.7 GeV, |? 7 (D*)| < 1.5} for the W range 
from 115 to 275 GeV is Acc = 11.4% for the above Monte Carlo program. 

We describe here the checks that are found to have a signihcant contribution to the sys¬ 
tematic error. For the acceptance, a systematic shift of 4-18% is estimated using different 
structure function parametrisations and 4-14% using HERWIG with MRSD'_/GRV HO as a 
different Monte Carlo program. Also, a systematic error of 8 % is found by varying the cuts on 
Pt{K, 7t) between 0.3 GeV and 0.7 GeV and on prijs) between 0.1 GeV and 0.25 GeV. Adding 
these errors in quadrature yields a total systematic error in the acceptance calculation of ^ 3 '^%. 
The systematic error on the number of events is 10%, estimated by using different background 
parametrisations to £t the signals. The systematic error of the luminosity measurement is 3.3%. 

Using the above quantities we measure an ep cross section for D*^ production, 
a{ep D*±A:) = a{ep D*+X) + a(ep ^ D*-X), of: 

a{ep —D*^X) = 32 ± 7{stat)~^‘^{syst) nb 

in the kinematic region {pt(D*) > 1.7 GeV, |? 7 (D*)| < 1.5} and 115 < W < 275 GeV. This 
cross section is valid for < 4 GeV^. The statistical error also includes the one due to the 
Monte Garlo statistics. 


In order to quote a cross section for charm production we need to correct for the fraction of 
events in which a charm quark pair fragments into D*"*" or D*“ as well as for the acceptance Aext 
of the kinematic region > 1.7 GeV, |? 7 (D*)| < 1.5}. The former is (52.0 ± 4.2)% 

and the latter is calculated by using PYTHIA with MRSD'_/GRV HO to be Ag^t = 13.7% for 
the region 115 < IF < 275 GeV. This extrapolation outside the kinematic region has a large 
uncertainty. In extrapolating pr(D*), the uncertainty is mainly due to the strong dependence 
on the rric value and for ? 7 (D*) it comes from the large differences between the different structure 
function parametrisations in the region |77(D*)| > 1.5. As a consequence, the systematic error 
of the product Acc ■ A^xt is very large. We have fixed rric to 1.5 GeV and quote the systematic 
error A{Acc ■ Aext) coming from the different structure functions and using HERWIG (SF and 
MG in Table respectively). Using a value of rric of 1.35 GeV (1.8 GeV) results in a shift of 
+25% (—40%) of the estimated cross section. 
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{W) 

(GeV) 

N 

Aec 

(%) 

-^ext 

(%) 

A(Acc 

SF 

MC 

a{ep ceX) 
(pb) 

Integrated 

$ 

aipip —> ccX) 
(ph) 

163 ± 16 

21 ± 7 

8.1 

16.2 


+54% 

0.23 ± 0.081^;^^ 

0.0367 

6.3 ± 2.2i;i;;) 

243 ± 24 

28 ±8 

22.4 

00 

bo 


+30% 

0.21 ± 0.061^:1^ 

0.0122 

TR9T5(21|^ 

198 ± 20 

48 ± 11 

11.4 

13.7 


+48% 

0.45±0.1ll)i:^^ 

0.0488 

9.1 + 2.21^;^ 


Table 1: Acceptances and cross sections 


We therefore estimate the ep charm prodnction cross section at ^/s = 296 GeV for < 
4 GeV^ in the range 115 < W < 275 GeV as: 

a{ep —*• ccX) = 0.45 ± 0.11 ^q;22 hb- 

This procednre was also carried ont dividing W into two ranges, 115 < W < 205 GeV and 
205 <W < 275 GeV. The (IV) for the events in these two ranges were 163 GeV and 243 GeV 
respectively. The error in W is dominated by the systematic nncertainty of the Jacqnet-Blondel 
method. The resnlts are shown in Table |l|. 


7.2 7 p Cross Section 


The 7 p cross section can be obtained from the corresponding ep cross section using the EPA 
formula: 

rymax rQmax o . o. . o. 

aep(s) = / dy ! dQ^ ■ Q^), 

•'ymin '^Qmin 


where 




a 1 
27r J/Q2 


[1 + (1 - vf 


2mly‘^ 


is the flux of transverse photons, Qmax — 4GeV^ and rUe is the electron mass. 

Since the median ^ 10“^ GeV^ is very small, we can assume the photons to be on-shell 
and therefore neglect the longitudinal contribution and the dependence of a^p. The yp 
cross section is then obtained as the ratio of the measured ep cross section and the photon 
flux factor 4) integrated over the and y range covered by the measurement. This procedure 
assumes that a^piW) is independent of y in the range of the measurement. As this dependence 
is not known a priori, the above calculation procedure was repeated assuming a proportional or 
logarithmic rise in W. An increase of 5% in the resulting cross section was found at (IV) = 163 
GeV and less than 2% at (IV) = 243 GeV. 


The estimated charm photoproduction cross section is (6.3 ± 2.2l|o) /ib at (fV) = 163 GeV 


we 


and (16.9 ± 5.2l8;’5 ) pb at (fV) = 243 GeV, assuming rric = 1.5 GeV (see Table |^). In Fig. || 
show our measurements for the total cc photoproduction cross section as a function of IV, in ad¬ 
dition to earlier measurements in hxed target experiments and NLO QGD calculations 
for rric = 1.5 GeV. The solid line represents the prediction with the MRSD(_/GRV HO struc¬ 
ture function parametrisation for proton/photon using = rric as the renormalisation scale. 
The shaded band represents the theoretical uncertainties coming from varying this scale in the 





























range 0.5 < HR/rric < 2. We also show the extreme predictions of MRSD'_/LAC1 (dashed line) 
and MRSDg/GRV HO (dash-dotted line) for = rric- We note that our measurement of the 
total charm photoproduction cross section at these high W values is in good agreement with 
those calculated using singular gluon distributions in the proton like MRSD'_ or GRV 
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8 Summary 


We have observed 48 D*(2010)^ mesons in the decay channel D° 7 r+ —(K“ 7 r+) 7 r+ (-|- c.c.) 

in photoproduction events from ep collisions at HERA. Gross sections have been determined 
for these events with < 4GeV^ and 115 <W < 275 GeV. The ep cross section for inclu¬ 
sive production is found to be (32 ± 7 ^ 7 ) nb in the kinematic region {pr(D*) > 1.7 GeV, 
| 77 (D*)| < 1.5}. Extrapolating outside this region and assuming a mass of the charm quark of 
1.5 GeV we estimate the ep charm cross sections to be a{ep —>• ccX) = (0.45 ± O.llt^^D hb at 
^/s = 296 GeV and 115 < IV < 275 GeV. The average yp charm cross section a{^p —> ccX) is 
found to be (6.3 ± 2.21 |q) ph at (IV) = 163 GeV and (16.9 ± 5 . 2 t 8 ^ 5 ®) ph at (PV) = 243 GeV. 
NLO QCD calculations using a gluon momentum density in the proton ~ at low x are 

in good agreement with the observed increase of the cross section by one order of magnitude 
when the jp centre-of-mass energy increases one order of magnitude with respect to previous 
low energy experiments. 


Acknowledgements 

We thank the DESY Directorate for their strong support and encouragement. The remark¬ 
able achievements of the HERA machine group were essential for the successful completion of 
this work, and are gratefully appreciated. We also gratefully acknowledge the support of the 
DESY computing and network services. We thank S. Frixione, M.L. Mangano, P. Nason, and 
G. Ridolh, for providing us with their results for the curves in Fig. 3. 


9 





References 


[1] A. Ali et al., in Proc. of the HERA Workshop, DESY (1987) 395. 

[2] R. van Woudenberg et al., in Proc. of the Workshop Physics at HERA, DESY (1991) 739. 

[3] E. Witten, Nucl. Phys. B120 (1977) 189; 

J.F. Owens, Phys. Rev. D21 (1980) 54. 

[4] M. Drees and F. Halzen, Phys. Rev. Lett. 61 (1988) 275; 

G.A. Schuler and T. Sjostrand, Phys. Lett. B300 (1993) 169. 

[5] L.E. Gordon and J.K. Storrow, Phys. Lett. B291 (1992) 320; 

D. Bddeker, Phys. Lett. B292 (1992) 164. 

[6] ZEUS Gollab., M. Derrick et ah, Phys. Lett. B297 (1992) 404. 

[7] ZEUS Gollab., M. Derrick et ah, Phys. Lett. B322 (1994) 287. 

[8] HI Gollab., T.Ahmed et al., Phys. Lett. B297 (1992) 205. 

[9] R.K. Ellis and P. Nason, Nucl. Phys. B312 (1989) 551. 

[10] S. Frixione et al., Phys. Lett. B308 (1993) 137. 

[11] P. Nason, S. Dawson and R.K. Ellis, Nucl. Phys. B303 (1988) 607; 

J. Smith and W.L. van Neerven, Nucl. Phys. B374 (1992) 36. 

[12] A.D. Martin, W.J. Stirling and R.G. Roberts, Phys. Lett. B306 (1993) 145. 

[13] M. Gliick, E. Reya and A. Vogt, Phys. Lett. B306 (1993) 391. 

[14] M. Drees and K. Grassie, Z. Phys. G28 (1985) 451. 

[15] P. Aurenche et ah, Z. Phys. G56 (1992) 589. 

[16] H. Abramowicz, K. Gharchula and A. Levy, Phys. Lett. B269 (1991) 458. 

[17] Particle Data Group, Review of Particle Properties, Phys. Rev. D50, (1994) 1173. 

[18] G. Alvisi et ah, Nucl. Instr. Meth. A305 (1991) 30. 

[19] N. Harnew et ah, Nucl. Instr. Meth. A279 (1989) 290; 

B. Foster et al., Nucl. Phys. B (Proc. Suppl.) 32 (1993) 181 and Nucl. Inst. Meth. A338 
(1994) 254. 

[20] A. Andresen et al. Nucl. Inst. Meth. A309 (1991) 101; 

A. Galdwell et ah, Nucl. Inst. Meth. A321 (1992) 352. 

A. Bernstein et ah, Nucl. Inst. Meth. A336 (1993) 23; 

[21] J. Andruszkow et al., DESY preprint, DESY 92-066 (1992). 


10 



[22] The ZEUS Detector, Status Report 1993, DESY 1993; 

W.H. Smith et ah, DESY preprint, DESY 94-183 (1994). 

[23] F. Jacquet and A. Blondel, Proc. of the Study for an ep Facility for Europe, ed. U. Amaldi, 
DESY 79/48 (1979) 391. 

[24] G. Marchesini et ah. Comp. Phys. Comm. 67 (1992) 465. 

[25] T. Sjostrand, Z. Phys. C42 (1989) 301 and in Proc. of the Workshop Physics at HERA, 
DESY (1991) 1405. 

[26] B. Andersson et ah, Phys. Rep. 97 (1983) 31. 

[27] T. Sjostrand, Comp. Phys. Comm. 39 (1986) 347; 

T. Sjostrand and M. Bengtsson, Comp. Phys. Comm. 43 (1987) 367. 

[28] C.F. Weizsacker, Z.Phys. 88 (1934) 612; 

E.J.Williams, Phys.Rev. 45 (1934) 729. 

[29] V.M. Budnev et ah, Phys. Rep. C15 (1975) 181. 

[30] W.K. Tung, Proceedings of International Workshop on Deep Inelastic Scattering and Re¬ 
lated Subjects, Eilat, Israel, 1994 World Sci., Singapore, to be published; 

H.L. Lai et ah, MSU-HEP-41024, Oct. 1994, to be published. 

[31] M.Gliick, E.Reya and A.Vogt, Z. Phys. C53 (1992) 127. 

[32] ZEUS Collab., M. Derrick et ah, DESY preprint, DESY-94-176. 

[33] OPAL Collab., G. Alexander et ah, Phys. Lett. B262 (1991) 341, updated in ICHEP94 
Ref. 0513, to be published in Z. Phys. C. 

[34] S. Frixione, M.L. Mangano, P. Nason and G. Ridolh, private communication. See also 
Nucl. Phys. B412 (1994) 225 and CERN-TH.7292/94, GEF-TH-4/1994, to be published 
in Nucl. Phys. B. 

[35] GIF Collab., M.S. Atiya et ah, Phys. Rev. Lett. 43 (1979) 414; 

BFP Collab., A.R. Clark et ah, Phys. Rev. Lett. 45 (1980) 682; 

SLAG HEP Collab., K. Abe et ah, Phys. Rev. D30 (1984) 1; 

EMC Collab., M. Arneodo et ah, Z. Phys. C35 (1987) 1; 

PEC Collab., M. Adamovich et ah, Phys. Lett. B187 (1987) 437; 

E691 Collab., J.C. Anjos et ah, Phys. Rev. Lett. 65 (1990) 2503; 

NA-14' Collab., M.P. Alvarez et ah, Z. Phys. C60 (1993) 53. 

[36] M. Chick, E. Reya and A. Vogt, Phys. Lett. B285 (1992) 285. 


11 




140 145 150 155 160 165 

AM (MeV) 


Figure 1: AM = M(K7r7rs') — M(K7r) distribution for DIS candidates : > 4GeV^ and 

Ve < 0.7. 
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Figure 2: a) AM distribution for photoproduction events having 1.80 < M{K7r) < 1.93 GeV: 
signal (dots) and wrong charge combinations (dashed histogram). The dashed line is a fit to the 
wrong charge background using the parametrisation A (AM — m^r)^, where 771^^ is the mass of 
the pion. The solid line is a fit to the distribution, parametrised as a sum of the same function 
for the background plus a Gaussian for the signal, b) {Kti) invariant mass distribution for 
those candidates with 142 < AM < 149 MeV. The htting function is the sum of a Gaussian 
and an exponential. 
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Figure 3: Total cc photoproduction cross section as a function of W. The solid dots are the 
ZEUS measurements and the open dots are earlier measurements from hxed target experiments. 
The inner error bars are the statistical and the outer ones the systematic errors. The solid line 
represents the central prediction of NLO calculations with MRSD'_/GRV HO parametrisations 
of the proton/photon structure function using /ir = rric (for rric = 1.5 GeV) as the renormal¬ 
isation scale. The shaded band represents the theoretical uncertainties coming from varying 
this scale in the range 0.5 < ^R/rric < 2. The dashed line represents the central prediction of 
MRSD'_/LAG1 and the dash-dotted line is MRSD[)/GRV HO. 
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